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Abstract 

This paper constructs a tree structure for the music rhythm using 
the L-system. It models the structure as an automata and derives 
its complexity. It also solves the complexity for the L-system. This 
complexity can resolve the similarity between trees. This complex- 
ity serves as a measure of psychological complexity for rhythms. It 
resolves the music complexity of various compositions including the 
Mozart effect K488. 
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1 Introduction 

This paper presents a computational model to interpret the 'Mozart effect' 
K488. This effect has been discussed extensively and seriously among both 
psychology and music perception societies using various experimental tech- 
niques [Rauscher et al. 1993]. Instead of experiments, this paper constructs 
a computation model to resolve the effect. This model starts with a tree 



structure for quantized rhythm beats based on the theory by Loguet-Higgins 
[Longuet-Higgins, 1987]. The tree is then modeled as an automata and its 
complexity is derived by way of the L-system [Prusinkiewicz and Linden- 
mayer 1990, Prusinkiewicz 1986]. The quantization tree will be briefly intro- 
duced in this section. The L-system will be introduced in the next section. 
The automata rewriting rule associated with the L-system will be included 
also in the next section. The similarity between trees by way of rewriting 
rules is defined in Section 3. The tree complexity is derived in Section 4. 
This complexity serves as a measure for the perception of musical rhythms 
[Desain and Windsor 2000; Yeston 1976] and resolves the effect. 

The key features of music perception and composition are rhythm, melody 
and harmony. Rhythm is formed through the alternation of long and short 
notes, or through repetition of strong and weak dynamics [Cooper and Meyer 
I960]. Because one metrical unit, such as a measure or a half note, can of- 
ten be divided into two or three sub-units (illustrated in Figures [H and ED , 
this rhythm is endowed with a clear hierarchical structure [Longuet-Higgins, 
1987; Lerdahl and Jackendoff 1983]. Note that Longuet Higgins grammar 
for rhythm is different from his musical parser for handeling performances, 
which is far more sophisticated than time-grid round-off. To represent the 
hierarchical characteristics of rhythm, we need to seek a system that pos- 
sesses such a nature. Fortunately, the plant kingdom is rich with branching 
structures, in which branches are derived from roots. In fact, the structure 
shown in Figure [T] is that of a binary tree. L-systems (Lindenmayer systems) 
[Prusinkiewicz and Lindenmayer 1990, Prusinkiewicz 1986] are designed to 
model plant development, see [McCormack 1993]. Therefore, it is natural 
to construct the rhythm representation by using an L-system [Prusinkiewicz 
1986; Worth and Stepney 2005]. A background of L-Systems applied to art 
and music is in the website in Reference. 

We will show how such a tree structure and its related parts can be 
constructed. We expect that the L-system can capture the rhythm nature. 
We now review the music tree by Loguet-Higgins. 

Figure [1] shows that in each level of a tree a half note is represented by a 
different metrical unit. In the highest level, the metrical unit is a half note; 
in the next level, the unit is a quarter note; in the lowest level, the unit is 
an eighth note. The total duration in each level is equal to a half note. This 
structure can be extended to a measure or more, a composition as in Figure 

m 

H. C. Longuet-Higgins introduced this kind of tree in the 1970s. The 
tree has been extensively studied both experimently and theoretically. Note 
that there are many other theories [Cooper and Meyer 1960; Yeston 1976; 
Lerdahl and Jackendoff 1983; Desain and Windsor 2000]. Since we prefer a 
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Figure 1: One metrical can be divided into two or three subunits. 



Figure 2: Musical tree of Beethoven's Piano Sonata No. 6, Mov. 3. 



computational approach to resolve the effect, we will not use their theories. 
In order to give computers the musicianship necessary to transcribe a melody 
into a score, he used tree structures to represent rhythmic groupings. In his 
theory of music perception, the essential task in perceiving the rhythmic 
structure of a melody is to identify the time of occurrence of each beat. 
Therefore, his theory can be applied to Western music with regular beats. In 
Western music, the most common subdivisions of each beat are into two or 
three shorter metrical units, and these shorter metrical units can be further 
subdivided into two or three units. Tracking from the start of a melody, when 
a beat or a fraction of a beat is interrupted by the onset of a note, it is divided 
into shorter metrical units. After this process of division, every note will find 
itself at the beginning of an uninterrupted metrical unit. The metrical units 
can be considered as the nodes of a tree in which each non-terminal node 
has two or three descendants. The terminal nodes for a beat are the shortest 
metrical units that the beat is divided into. Every terminal node in the tree 
will eventually be attached either to a rest or to a note sounded or tied. It 
is natural to include and elaborate rests in the tree model as those done in 
[Longuet-Higgins, 1987]. 

We will employ the perception factors discussed by Longuet-Higgins, such 
as tolerance, syncopation, rhythmic ambiguity, regular passages, to construct 
L-systems for rhythms. 
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2 Rhythm represented with tree structure 



A rhythmic tree as described above is a tree of which each subtree is also 
a rhythmic tree. Each tree node has two or three children (branches or 
descendants). Each node in the tree represents the total beat duration that 
is equal to the sum of all those of its descendants. The root node has a 
duration length that is equal to the length of the whole note sequence. 

Note that when we attempt to split a note sequence into two subsequences 
with equal duration lengths, we usually obtain two unequal length subse- 
quences. This is because a note connecting the two subsequences has been 
split into two submetrical units. The preceding portion belongs to the pre- 
ceding subsequences and the later portion belongs to the later subsequence. 
We will mark those units to identify their subsequences. 

These two subsequences represent two different subtrees of the root node. 
We further divide each subsequence into sub-subsequences, which are also 
rhythmic trees, and so on. This dividing process is completed when a tree 
node contains a single note. This single note may possibly be the one which 
has been split into two portions. This is in some sense similar to an algorithm 
for note quantization and is a standard practice in MIDI rendering of music. 
In practice we will quantize notes using the finest note among dotted notes 
(e.g., 1/4, 1/8, 1/16, dotted notes, etc.) without loosing most of the interest- 
ing details. We plot two such trees in Figures [2] and El The notes shown in 
Figure [3] are part of the whole tree of the beginning of Rachnaminoff's Piano 
Concerto No. 3, Movement 1. The two notes in the rectangle have been split 
using our rhythmic tree process. 



Figure 3: Musical tree of Rachmaninoff's Piano Concerto No. 3, Mov. 1. 
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2.1 Rhythm represented with L-system 

To express the hierarchical characteristics of rhythm, we need a data struc- 
ture that possesses such a hierarchical nature. Fortunately, the plant king- 
dom is dominated by branching structures, in which branches are derived 
from roots. L-systems are designed to model plant development. Therefore, 
it is practicable to construct a rhythmic representation by using L-systems. 

The Lindenmayer system, or L-system for short, was introduced by the 
biologist Aristid Lindenmayer in 1968 [Lindenmayer 1968]. It was conceived 
as a mathematical theory of plant development. The central concept of the 
L-system is rewriting. In general, rewriting is a technique used to define 
complex objects by successively replacing parts of a simple initial object, 
using a set of rewriting rules or productions. 

The L-system is a new type of string-rewriting mechanism. The essential 
difference between Chomsky grammars and L-systems lies in the technique 
used to apply productions. In Chomsky grammars, productions are applied 
sequentially, whereas in L-systems, they are applied in parallel and simulta- 
neously to replace all the letters in a given word [McCormack 1993]. 

This difference reflects the biological motivation of L-systems. Produc- 
tions are intended to capture cell divisions in multi-cellular organisms, where 
many divisions may occur at the same time. Moreover, there are languages 
which can be generated by context-free L-systems but not by context-free 
Chomsky grammars. 

Here we introduce the turtle graphical interpretation of L-systems. Sup- 
pose that there is a turtle crawling on a plane. The state of the turtle is 
defined as a triplet (x, y, a), where the Cartesian coordinates (x, y) represent 
the turtle's position, and the angle a, called the heading, is interpreted as 
the direction in which the turtle is facing. Given the step size d and the 
angle increment 5, the turtle can respond to commands represented by the 
following symbols: 

F Move forward a step of length d. The state of the turtle changes 
to (x new , y new , a) , where x new = x + dcosa and y new = y + dsina. 
Draw a line segment between points (x,y) and (x new ,y new ). 

f Move forward a step of length d without drawing a line. 

+ Turn left (counterclockwise) by angle S. The next state of the turtle 
is (x,y,a + 5). The positive orientation of angles is counterclock- 
wise. 

Turn right (clockwise) by angle 5. The next state of the turtle is 
(x,y,a-8). 
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For a string composed by the above symbols, the turtle will crawl accord- 
ing to commands indicated in a given order. The turtle interpretation of this 
string is the figure drawn by the turtle. To draw branches of a tree, we need 
two more symbols: 

[ Push the current state of the turtle onto a pushdown stack. The 
information saved on the stack contains the turtle's position and 
orientation, and possibly other attributes such as the color and 
width of the lines being drawn. 

j Pop a state from the stack and make it the current state of the 
turtle. No line is drawn although, in general, the position of the 
turtle changes. 

These two symbols enable the turtle to go back to the root after drawing 
a branch so that it can draw other branches originating from the same root. 
This kind of L-system is called the bracketed L-system, and we call strings 
that represent trees bracketed strings. 

2.2 Rewriting rules for rhythmic trees 

In music, a longer metrical unit can be replaced by the combination of several 
shorter metrical units, such as a bar filled with several notes. Given a rhythm, 
the direct way to represent it with an L-system is to construct rewriting rules 
that replace longer metrical units with shorter ones. Take the dotted half 
note for example; it can be replaced by 3 quarter notes or by the combination 
of a quarter note, a dotted quarter note and an eighth note, of a half note 
and a quarter note, and so on, see Figure HI In Figure HI 4 rewriting rules are 
shown that can be used to rewrite a longer metrical unit (a dotted half note) 
into several shorter metrical units (a quarter note, dotted quarter note, and 
so on). 



Figure 4: Different combinations of a dotted half note. 
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Thus, we can regard the whole score as the longest metrical unit and con- 
struct rewriting rules for it, recursively from the longest unit to the shortest 
metrical unit. These rewriting rules can represent the rhythm because we 
can regenerate the whole score by using these rules. We will not elaborate 
on the rewriting rules for a score [Lee and Liou 2003]. The rewriting rule 
system is formally equivalent to the bracketed L-system. 

Similar to a score, a rhythmic tree can also be represented by means of 
rewriting rules. After we build a rhythmic tree from a score by using the 
technique described in Section 2, we may apply rewriting rules to each node 
in the rhythmic tree. As a half node can be rewritten as 2 quarter nodes, 
each node can be rewritten as all the subtrees attached to it. The rewriting 
rules for our tree nodes always generate 2 subtrees because the technique 
described in section 2 always generates 2 subtrees (children) for each node. 
Note that similar qualitative results of this work can be obtained with 3 
or more subtrees. We will focus on 2 generated subtrees in this work. In 
Figure [5j we represent a subtree (or a metrical unit) using the rewriting rule 
P — > LR, where P denotes the subtree we want to represent using a set of 
rewriting rules, L denotes its left subtree, and R denotes its right subtree. 
We call L and R the nonterminals. With this schema, we can write the 
rewriting rule for the tree shown below. 



Figure 5: Using rewriting rules to represent a rhythmic tree. 

In Figure El Ll denotes the left subtree's left subtree, and Rl denotes the 
right subtree's left subtree. R Rl and R Rr are similar. Thus, the rewriting 
rules for the tree shown in Figure E] are P — ► L l L r R l R Rl R Rr . 

2.3 Bracketed strings for a rhythmic tree 

As discussed above, bracketed strings can represent a hierarchical structure, 
such as an axial tree. Thus, bracketed strings may provide a suitable data 
structure for representing rhythm. We know that a half note can be divided 
into two quarter notes, and this fact can be represented by a tree structure, 
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Figure 6: Bracketed strings for two trees. 



which has a parent node and two child nodes, as shown in Figure |6] (a). 
The bracketed string of the half note and its binary branches in Figure El (a) 
are F[+F][-F]: the first F is the command for tracing the root; [+F] is the 
command for tracing the left branch, and [-F] is the command for tracing 
the right branch. In Figure [6] (b), there is another example for a dotted half 
note. 



Figure 7: Bracketed string for Beethoven's Piano Sonata No 6, Mov. 3. 

Since we focus here only on rhythmic trees, we can simplify the brack- 
eted string representations. First, our rhythmic trees have only 2 subtrees. 
Second, the 'F' notation for a rhythmic tree is trivial. With these two char- 
acteristics, we may omit the 'F' notation from the bracketed string and use 
only four symbols, {[, ], -, +}, to represent rhythmic trees. In our cases, 
'[...]' denotes a rhythmic (sub)tree where '. . . ' indicates all the bracketed 
strings of its subtrees. '-' indicated that the next '[•••]' notation for a tree is 
a left subtree of the current (sub)tree, and '+' indicates that the next '[...]' 
notation is a right subtree. In Figure [H we list the simplified rules for the 
subtrees shown in Figure [71 

Note that this rhythmic tree has no middle subtree in each node. In our 
model, we always choose to divide a metrical unit into exactly two shorter 
metrical subunits. We will assume that all the rhythmic trees we generate 
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Figure 8: Bracketed strings for each node of the rhythmic tree shown in 
Figure [7J 

are such binary trees. To simplify the snalysis, we also assume that the notes 
are monophonic. Another example in Figure [9] is a bracketed string for the 
tree of Rachmaninoff's Piano Concerto, as shown in Figure [31 Here we list 
all symbols in the bracketed string. 



Figure 9: Bracketed string for Rachmaninoff's Piano Concerto No. 3, Mov.l. 



2.4 Rewriting rules for bracketed strings 

Now, we know how to use rewriting rules to represent a tree, and we also 
know how to represent a tree with a bracketed string. We can also use 
rewriting rules to generate bracketed strings. In rewriting rules for rhythmic 
trees, we write P — > LR for a tree having left and right subtrees. Note that 
we call L and R the nonterminals. Such a tree will have a bracketed string as 
follows: [[-F. . . ][+F. . .]]. It is clear that '[-F. . .]' represents the left subtree, 
and that '[+F. ..]' represents the right subtree. Therefore, we can replace 
the rewriting rules with 

P -> [-FL][+FR] 
L - 
R -> 

where '. . . ' is the rewriting rule for the bracketed string of each subtree. 
In this way, we do not have to write L for a left subtree and R for a right 
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subtree; the orientation is already described in the bracketed string £ -F' and 
£ +F'. Thus, we do not have to write words such as 'Rr l \ 1 Rr r \ etc. Of 
course, we may still use such recursive subscript representations for rules for 
the sake of readability. In Figure [TO], we show the rewriting rules for the 
bracketed string of the tree in Figure El 



Figure 10: Rewriting rules for the bracketed string of a rhythmic tree. 



There are "nulls" in the rules. We use "null" to represent a terminal 
(or a tree node that doesn't have any child subtree). For such null-subtree 
rewriting rules we simply ignore the nulls. The new rewriting rules without 
trivial nulls are as follows: 

P -> [— FTl][+FTr] 
T L - [-F][+F] 
Tr - [-F][+FT RR ] ■ 
Trr - [-F][+F] 
Note that there are two identical rules in the above rewriting rules: Tl and 
Trr. This redundancy raises a question: Can we combine them to simplify 
these rules? Doing so will not harm the whole structure if the redundant 
rules contain only null subtrees. We will show in the following what will 
happen if the rules do not contain only null subtrees. Assume that we have 
the following rules: 



p 


> [-FT L ][+FT R 


T L 


► [-F][+F] 


Tr 


► [-F][+FT RR ] 


Trr — 


► [-F][+FTrrr 


Trrr — 


► [-F][+F] 



These rules can generate exactly one bracketed string and, thus, exactly 
one rhythmic tree. All these rules form a rule set, which represents a unique 
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rhythmic tree. It is clear that T R and T RR are almost the same. The only 
difference is that one of the subtrees is T RRR , and that the other is T RR . But 
they have the same structure: they both have a right subtree and do not 
have a left subtree. We can use this similarity to explore the characteristics 
of a music work, and from a composer's works, we can explore his or her 
characteristics. We will define two terms to express the similarity between 
two rewriting rules. 

3 Homomorphism and isomorphism of rewrit- 
ing rules 

We will now study some characteristics of rewriting rules to extract similar 
sections of trees. We first define the similarity between two sections as follows: 

DEFINE : Homomorphism in rewriting rules. 

Rewriting rule R\ and rewriting rule R2 are homomorphic if and only 
if they have the same structure. Their corresponding rhythmic trees both 
have subtrees in corresponding positions or both not. That is, ignoring all 
nonterminals, if rule R\ and rule R2 generate the same bracketed string, then 
they are homomorphic by definition. 

DEFINE : Isomorphism on level X in rewriting rules. 

Rewriting rule R\ and rewriting rule R2 are isomorphic on depth X if 
they are homomorphic and their nonterminals are relatively isomorphic on 
depth X — 1. Isomorphic on level indicates homomorphism. 



Figure 11: Rhythmic tree. 

Here, we will use an artificial example to clarify these definitions. In 
Figure [HI we name the tree rooted at A, B, C, and D, respectively, tree A, 
tree B, tree C, and tree D. Tree A is homomorphic to tree B and tree C, but 
tree A is not isomorphic to tree D. Tree A is isomorphic to tree C on depth 
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2, but they are not isomorphic on depth 3. Tree B is isomorphic to tree C on 
depth and 1, but not on depth 2. D is not isomorphic to any other trees, 
nor is it homomorphic to any other trees. 

By using bracketed strings, we can obtain a much clearer definition of 
homomorphism. All homomorphic rewriting rules generate the same brack- 
eted strings when all nonterminals are ignored. Using the following rewriting 
rules, we find that P and T L are homomorphic because if T L and T R are ig- 
nored, they generate exactly the same bracketed string, [-F][+F]. We also find 
that P and Trr are homomorphic because if Trrr, Tl and Tr are ignored, 
they generate the same bracketed string, [-F] [+F] : 



In fact, in this example, all five rewriting rules are homomorphic to each 
other. But if we add a sixth rule, Trrrr —>■ [— F], then it will not be homo- 
morphic to any of the other rules. 

Once we define the similarity between rules, we can classify all the rules 
in the set into different subsets based on their similarity. Rules that belong 
to a class are all isomorphic to each other on depth X. All the rules' names 
are replaced with the names of the classes to which the rewriting rules be- 
long. After such name conversion, every new rewriting rule represents more 
rhythmic trees than it was. These new rewriting rules set can now generate 
more rhythmic trees, including the original one. Note that the definition 
may include certain over-generalization cases. We show one case in the last 
Section. 

We know that these rules can generate the original bracketed string, but 
it can actually do more than that. In fact, it can generate an infinite number 
of bracketed strings. After performing classification, we obtain not only 
a new rewriting rule set but also a context free grammar, which can be 
converted into an automata. In this case, we say that this rhythmic tree can 
be converted into an automata, which can generate many other rhythmic 
trees that have similar characteristics. We list the rewriting rules in Table 
1 for the previous example shown in Figures [3] and [91 Rule Other —>■ null 
is ignored, there are 22 such rules in the tree we list only one for simplicity. 
The classification of the rules is listed in Table 2. 




[-FT L ][+FT R ] 

[-F][+F] 
[-F][+FT RR ] 

[-F][+FT RRR ] 
[-F][+F] 
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P — LR 


R — > R L R R 


L — L L L R 


Rl^ RllRlr 


LllLlr 


Rll^ Rlll 


+F] 




Lll^ [-F] 


[+F] 


Rlll— > [— F 


[+F] 




Llr^ Llrl[+F] 


Rlr^ [-F][+F] 


Llrl^ [-F][+F] 


Rr^ RrlRrr 


Lr^ LrlLrr 


Rrl^ Rrll 


[+F 




Lrl^ Lrll[+F] 


Rrll- [-F][+F] 


Lrll^ [-F][+F] 


Rrr^ [~F]R rrr 


Lrr^ [— F 




Rrrr^ [-F][+F] 


Lrrr^ [— F][+F] 


Other — > null 



Table 1. Rewriting rules for the rhythmic tree in Figure 3. 



Classification 
of rules 


Isomorphic 
Depth #0 


Isomorphic 
Depth #1 


Isomorphic 
Depth #2 


Isomorzphic 
Depth #3 


Class #1 


(21)C 1 -> C 2 C 2 


(3)0!- Cid 
(2)C 1 -C 2 C 8 

(1)0^ 0402 


(1)C 1 -C 2 C 4 

(l)C 1 -C 8 C 4 


(^d-CaCa 


Class #2 


terminal 
(22) C 2 ^ null 


(4)C 2 ^ C 4 C 5 


(1)C 2 - C 7 C 5 


(1)C 2 - C 8 C 4 


Class #3 




(2)C 3 ^ C5C4 


(1)C 3 ^C 5 C 7 


(l)Cg- C9C4 


Class #4 




(8)C 4 — C5C5 


(2)C 4 ^ C 5 C 6 


(2)C 4 ^ C 5 C 6 


Class #5 




terminal 
(22)C 5 ^ null 


(4)C 5 ^ C 7 C 8 


(4)C 5 ^ C 7 C 10 


Class #6 






(2)C 6 ^ C 8 C 7 


(2)C 6 ^ CioC 7 


Class #7 






(8)C 7 ^ C 8 C 8 


(8)C 7 — C10C10 


Class #8 






terminal 
(22)C 8 ^ null 


(1)C 8 - C 7 C 5 


Class #9 








(1)C 9 - c 5 c 7 


Class #10 








terminal 
(22)C 10 ^ null 



Table 2: Classifying based on the similarity of rewriting rules. 



In Table 1, rules such as Rrrr — [— F][+F] , L RRR — [— F][+F] are 
assigned to Class 2. There are eight such rules before classification, so we 
write '(8)6*2 — [— F][+F]'. Similar rules such as Lr — LrlLrr, P — LR, 
R — > RlRr are isomorphic on depth 0, and there are 11 such rules. They are 
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assigned to Class 1. Class 3 and Class 4 are obtained by following a similar 
classification procedure. Note that this section also presents a new way to 
convert a context-sensitive grammar to a context-free one. 

4 Rhythmic complexity 

After we list the rewriting rules for a rhythmic tree and classify all those rules, 
we attempt to explore the redundancy in the tree (the hidden structure in the 
beats) that will be the base for building the cognitive map [Barlow 1989]. To 
accomplish this, we compute the complexity of the tree which those classified 
rules represent. We know that a classified rewriting rule set is also a context 
free grammar, so we can define the complexity of a rewriting rule set as 
follows: 

DEFINE : Topological entropy of a context free grammar. 

The topological entropy K of a CFG (Context Free Grammar) can be 
evaluated by means of the following three procedure [Kuich 1970; Badii and 
Politi 1997]: 

(1) For each variable Vi with productions (in Greibach form), 

Vi — > UiUu, t i2 Ui2, UkJJiki, 
where *i2, *i3— *ifcj are terminals and {U a , U i2 , ...U iki } are 
non-terminals. The formal algebraic expression for each variable is 

Vi = Y2j = i UjUij. 



(2) By replacing every terminal Uj with an auxiliary variable z, one ob- 
tains the generating function 

where iVj(n) is the number of words of length n descending from V$. 

(3) Let N(n) be the largest one of iVj(ra), N(n) = max{N i (n) ,over all i}. 
The above summation series converges when z < R = e~ K ". The topological 
entropy is given by the radius of convergence R as 

K = -1n R. 
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However, we have found that this definition is slightly inconvenient for 
our binary tree case. Thus, we rewrite it as follows: 

DEFINE : Generating function of a context free grammar. 

Assume that there are n classes of rules and that each class Cj contains 
m rules. Let Vj G {Ci, C 2 , C 3 , C n }, Uy G {Rij,i = 1~ nj = 1~ nj, and 
a.ijk G {x : x = 1~ n}, where each has the following form: 

— > Kin Kiia 

^£2 ~ ► K 4 2lKi 22 

The generating function of Vi, Vi(z) , has a new form as follows: 

Z^<j=l "s<J 

If Vi doesn't have any non-terminal, we set Vi(z) = 1. With this function, 
we can define the complexity of the rhythmic tree below. 

DEFINE : Complexity of rhythmic tree [6]. 

After formulate the generating function Vi(z), we intend to find the largest 
value of z, z max , at which Vi(z max ) converges. Note that we use V\ to denote 
the rule for the root node of the rhythmic tree. After obtaining the largest 
value, z max , of Vi(z), we set R = z max , the radius of convergence of V\(z). 
We define the complexity of the rhythmic tree as Kq = — In R. 

We use the simple example in Tables 1 and 2 (or Figure [3]) to show the 
computation procedure of the complexity. According to our definition the 
given values for the class parameters are {n = 5, ri\ = 4, n 2 = 1, ra 3 = 1, 714 = 

1, ra 5 = l,nu = 3,n 12 = 2,n 13 = 1, n 1A = l,n 21 = 4, n 31 = 2,n 41 = 8,n 51 = 
22, dm = l,an2 = 1,0121 — 2, 0122 = 3,ai3i = 2, 0132 = 4, aui = 4, au 2 = 

2, ct 2 ii = 4, 0212 = 5, 0311 = 5 , 0312 = 4, 0411 = 5, 0412 = 5, a^u = 2, 0512 = 3}. 

Substituting these values in the equation, we have V 5 (z') = 1 and V±{z') = 
z' directly. Then we obtain the formulas for Vs(z), V2(z), and Vi(z) succes- 
sively. They are 



V 3 (z> 
V 2 (z> 
VAz' 



(Epii n 3p z'Vg 3pl (z')Vg 3p2 jz')) _ z 'x(2xV 5 (z')xV 4 (z')) _ , 2 
(Ep£i n 2 pz'Vg 2pl (z')Va 2p2 (z')) _ z 'x(4xV 4 (z')xV 5 (z')) _ , 2 



{Yr v l-yn Xv z>V aivl {z')V aXv2 {z>)) 

9=1 n M 

z'x(3xV 1 (z') 2 +2xV2(z')xV3(z')+lxV2(z')xV i (z')+lxV4(z')xV 2 (z')) 

7 

3z' x Vx (z') 2 +2 x (z'f+2 x (z') 4 
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Rearranging the above equation for Vi(z), we obtain a quadratic equation 
for V^z') 

Kz'mz') 2 - V^z') + f {{z'f + (z') 4 ) = 
Solving Vi(z'), we obtain the formula 

v(yl \ _ i±yi-i((^) 5 +(^) 4 ) 

V ^ Z ) ~ Wj/7 • 

The radius of convergence, R, and complexity, K = — In R, can be ob- 
tained from this formula. 



5 Implementation and Examples 

In order to compute the complexity of a rhythmic tree, we have to determine 
R, the radius of convergence of the rhythmic tree's rewriting rule set. We 
devise strategy to judge whether the function Vi(z') is convergent or divergent 
for a given value of z'. We construct an iteration technique to compute the 
value of this generating function. To facilitate the computation, we rewrite 
the generating function as follows: 

Epli n ip z'V^- 1 1 (z')V^~ 2 1 (z') 

1. 

Here we use superscript m in the variable V™(z') to represent the iteration 
count. Starting with V^(z') in each iteration, we calculate a new value, 
Vl(z'). Then we calculate V?(z'), V?(z'), ... , and V™{z') successively, 
where m is some positive integer number. When VJ n ~ 1 (z') is equal to V j m (z') 
for all rules, this means that V™'{z') cannot be improved anymore, we reach 
convergence. Therefore, z', the number we want to judge, is not the radius of 
convergence for the rules set but is smaller than the radius of convergence. In 
our simulations, we set m = 1000. This means that if V i m (z') is not divergent 
for m < 1000, then we judge z' to be convergent. 

Once we can judge whether Vi(z') is convergent or divergent at a number 
z', we can test every real number between and 1 to find the number that is 
right on the border of the convergent region and use this number to calculate 
the radius of convergence. We may apply some advanced techniques to search 
for the radius of convergence, such as binary searching between and 1. This 
is exactly the technique we use in our algorithm. 



vr(z') = 

KV) = 
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Now we will present a practical example. We use Beethoven's Piano 
Sonatas Nos. 1 to 32 and Mozart's Piano Sonatas Nos. 1 to 19 as an exam- 
ple, and show their complexity. We list the complexity of each piano sonata 
by Mozart in Figures fT2lfT3l In these figures, we use two different isomorphic 
depths, 1 and 3, to compute the complexity. From the figures we can see 
that the complexity is high for both composers. When we use higher depth 
isomorphism to classify rules, the complexity will decrease. This is because 
when we use higher depth isomorphism, redundancy between rules will de- 
crease so the complexity will also decrease. Eventually the complexity will 
decrease to zero for the highest depth isomorphism. Conversely, lower depth 
isomorphism brings more rules in a class; redundancy between rules will in- 
crease and the number of classes will decrease. If the depth of isomorphism 
is too low, the rules set will become too simple, thus the complexity will also 
become lower. We may compute the complexity for different depths to see 
the differences. 

Beethoven and Mozart's work have similar complexity, but Beethoven's 
is slightly higher than Mozart's. Both their complexity isomorphic on level 
2 are the highest. When we use the high level isomorphism to classify rules, 
the complexity of rules will decrease. Reversely, the low level isomorphism 
collects many rules in a class; redundancy between rules will increase and 
the number of classes will decrease. If the level of isomorphism is too low, 
the rules set will become too simple, thus the complexity will also become 
lower. We can try different level to see its complexity, and pick up the level 
with highest complexity. 

We tested a well-known music work studied by Rauscher et al. [1993]. 
Almost all the previous studies on the Mozart Effect have focused on a single 
piece of music, the Sonata for Two Pianos in D Major (K448). We have 
computed its complexity and found that it is generally higher than that of 
other sonatas by Mozart, see Figures fT2lfT3l 

6 Discussion 

We have constructed the complexity for the L-system. This complexity re- 
sembles, in some sense, the redundancy [Pollack 1990; Large et al. 1995; 
Chalmers 1990]. This complexity can facilitate many other studies such as 
bio-morphology, DNA analysis, gene analysis and tree similarity. 

We closely followed the ideas of Barlow [Barlow 1989] and Feldman [Feld- 
man 2000] to design this model. In his work, Barlow wrote that: 'Words 
are to the elements of our sensations like logical functions to the variables 
that compose them. We cannot of course suppose that an animal can form 
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an association with any arbitrary logical function of its sensory messages, 
but they have capacities that tend in that direction, and it is these capaci- 
ties that the kind of representative schemes considered here might be able to 
mimic." Human perception sometimes bases on external world's information 
redundancy If we can extract any rules or patterns from a certain object as 
part of our cognition map for that object, it will be easy to memorize or com- 
prehend it. In our model, rhythms resemble the words; trees resemble the 
logical functions; classes resembles rules and patterns; complexity resembles 
redundancy. 

Man is not inherently musical, the distinguished scientist Newton claimed; 
natural singing is the sole property of birds. In contrast to our feathered 

friends, humans perform and understand only what they taught This 

is why humans listen to music by training. One needs such redundancy to 
comprehend the music words. 

But how can we pinpoint the rules or patterns in a music work, or even 
in a simple rhythm that may be formless? As an attempt, we have defined 
homomorphism and isomorphism so as to characterize the similarity between 
sections of different rhythmic trees. But there still exist questions about the 
psychological implications of these characteristics, such as the depth of iso- 
morphism. The proposed model can enable us to measure the psychological 
complexity [Feldman 2000] of rhythms. In our studies, we have found that 
different depths of isomorphism produce varying degree of complexity. If a 
rhythm is very simple, its complexity will be 0. The same situation also oc- 
curred when we used isomorphism with a very high depth value to compute 
the complexity of Mozart's and Beethoven's piano sonatas. In general, the 
results confirm our intuition about these musical rhythms. 

Note that not all music can be properly approached with binary struc- 
tures. If the structure of a music piece is ternary, we expect that the com- 
puted complexity will be higher than it is in reality. 

We define the similarity between tree structures in Section 3. Finding 
similarity between rules and classifying them in different subsets are in some 
sense similar to fractal compression, see the website in Reference. This could 
be an alternate way to configure rhythmic complexity. We are still work- 
ing on this. We are also working on an extension of the model to incorpo- 
rate the rhythmic complexity for polyphonic music, superposition of different 
rhythms, tempo variation, grace notes, supra and irregular subdivisions of 
the beat (e.g. triplets, quintuplets,...). 
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Figure 12: Mozart's 19 Piano Sonatas, using isomorphic depth 1. 



Figure 13: Mozart's 19 Piano Sonatas, using isomorphic depth 3. 
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